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Abstract 

An non-autonomous system is proposed to model the seasonal pat- 
tern of dengue fever. 

We found that an approximate threshold condition for infection 
persistence describes all possible behavior of the system. 

As far as we know, the kind of analysis here proposed is entirely 
new. No precise mathematical theorems are demonstrated but we give 
enough numerical evidence to support the conclusions. 
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1 Introduction 

In a previous paper ilj, in which we tried to understand the phenomenon 
of dengue overwintering, we discovered an interesting threshold condition 
that allows the complete qualitative understanding of the behavior of a non- 
autonomous system. 

To motivate the reader we briefly describe the phenomenon we studied in 

in 

In subtropical regions dengue fever, a mosquito transmitted disease, shows 
a resurgent pattern with yearly epidemics, which starts typically in the 
months characterized by heavy rains and heat, peaking some three or four 
months after the beginning of the rainy season. In the dry months the num- 
ber of cases drop essentially to zero due to the virtual disappearance of the 
vector. Since the infection reappears for some years in the same regions, it 
is natural to ask how the virus survives the dry season. 

In order to model those seasonal patterns of the disease, we proposed a 
non-autonomous system, described below. 

The model describes the dynamic of dengue in its three components of 
transmission, namely, human hosts, mosquitoes and their eggs (the latter 
includes the intermediate stages, like larvae and pupae). These populations, 
in turn, are divided into susceptible humans, denoted Sh, infected humans, 
Ih, recovered (and immune) humans, Rh, total humans, = Sh+Ih+Rh, 
susceptible mosquitoes, Sm, infected and latent mosquitoes, Lm, infected and 
infectious mosquitoes, Im, non-infected eggs, Se, and infected eggs, Ie- 

The model's dynamics is described by the set of equations. 
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^ = -ablM^-I^HSH + rnNHil-^) 

^ = ablM^-{l^H + »H + lH)lH 

^ = P5 (C5 -(^5 sin (27r/i + 0)) -5^^(05-^5 sin (27r/i + 0)) 

—fJ'M^M — O-Smj^ 

P7 (c/ - rf/ sin {271 ft + (j))) IeO (c/ - dj sin (27r/t + 0)) 

^ = [rM5M+(l-5)rM/M](l-^^)- 

HeSe - Ps {cs - ds sin {27: ft + tt)) SeO {ce - dE sin {27: ft + tt)) 



If = 9rMlM{l-^-^)-f^ElE' 

Pi {ci - di sin {27: ft + 0)) /g^ (c/ - sin (27r/t + 0)) 

(1) 

Let us briefly describe some features of the model. 

We begin by describing the first three equations of the model. 

Susceptible humans grow at the rate rj:fiV//(l — ^) — (Ij^Sh, where rn 
is the birth rate, fi^^ is the natural mortality and k.H is the human carrying 
capacity. Note that we are assuming that close to the carrying capacity 
the human population growth is checked by a reduction in the birth rate. 
Alternatively the control of the population could be done by assuming an 
increase in the mortality rate, but the net result would be qualitatively the 
same. Those susceptible humans who acquire the infection do so at the rate 
ahlMjf^, where a is the average daily biting rates of mosquitoes and 6 is 
the fraction of infective bites inflicted by infectious mosquitoes Im- Infected 
humans, Ir may either recover, with rate 7, or die from the disease, with 
rate {^fj + an)- Recovered humans remain so for the rest of their lives. 

The fourth, fifth and sixth equations represent the susceptible, latent and 
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infectious mosquitoes sub-populations, respectively. Susceptible mosquitoes 
varies in size with a time-dependent rate 

Ps {cs - ds sin (27r/i + 0)) SeO {cs - dg sin (27r/i + (fy) - HmSm- 

The term jij^ is the natural mortality rate of mosquitoes. The term psSe is 
the fraction of eggs present at time t, and which survived the development 
through the intermediate stages (larvas and pupas). The time-dependent 
rate (q — di sin (27r/t + (t>)) {ci — di sin (27r/t + (p)) {i = S, I) simulates the 
seasonal variation in mosquitoes production from eggs, assumed different for 
infected and susceptible eggs, for generality. By varying Cj and di , {i — S, I), 
we can simulate the duration and severity of the winters (/ — 1/365 days~^ 
and so it fixes one cycle per year). The Heaviside ^^-function (a step function 
that is equal to zero when the argument is less than zero and one when the 
argument is greater or equal to zero) 9 (q — sin (27r/t + 0)) prevents the 
term 

{ci - di sin (27r/t + 0)) ^ (q - di sin {2% ft + 0)) {i = S, I) 
from becoming negative. If q is smaller than di, then the winter is long 
and severe. On the other hand, if q is greater than di, then the winter is short 
and mild. Susceptible mosquitoes who acquire the infection do so at the rate 
o-Smj^, where a is the average daily biting rates of mosquitoes, and became 
latent. A fraction of the latent mosquitoes survives the extrinsic incubation 
period with probability e'~^^''^' and become infectious. Therefore, the rate 
of mosquitos becoming infectious per unit time is e~^^'^' aSuit ~ tj) ^"^^'^ ■ 
The term 



represent vertical transmission, that is, the rate by which infected eggs be- 
come infectious adults. Infected mosquitoes die at the same rate jij^ as the 
susceptible ones. 

The seventh and eighth equations represent the dynamics of susceptible 

and infected eggs, respectively. 

In the seventh equation, the term 



represent the birth rate of susceptible eggs born from susceptible mosquitoes 
with rate 



Pi {ci - dj sin (27r/t + </>)) IeO {cj - di sin (27r/t + ({))) 





4 



and from a fraction {1 — g) of infected mosquitoes, with rate 



(1 



The term tm 



1 — ^ ^, ) represents the density-dependent rate of eggs 




birth rate. Once again we choose a density dependence on birth rather than 
on death. Ahernatively the control of the population could be done by 
assuming an increase in the mortality rate fi^, but the net result would be 
qualitatively the same. Finally, in the last equation the term 



as already mentioned, is the fraction of the hatched infected eggs which 
evolves to infectious adults. 

2 An approximated threshold condition 

In the first part of this section we deduce a threshold condition for epidemic. 
The intuition behind the procedures is discussed later on. 

In order to deduce the threshold condition for epidemic we replace the 
non-autonomous system ([Q) by a autonomous one, by regarding the time on 
the right side of the system ([Q) as a parameter and then carry out a local 
stability analysis. We linearize the second, the fifth, the sixth and eighth 
equations of the autonomous system around a small amount of disease in, 




represents the rate by which infected eggs 



TOW and the term 



Pj {cj — dj sin {271 ft + 0)) IeO (c/ — sin {2tt ft + 0)) , 
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Im, iM and tE- 



dt N 



Pi {cj — dj sin ($)) (c/ — c// sin ($)) 

If = ^^-M (l - ^) - /X^Zi.- 

(c/ — (i/ sin ($)) isO (c/ — sin ($)) 

where $ = 27r/t + (p. 

We then examine the stabihty of the trivial solution of system (j2I), that 
is, z^; = 0, Im = 0, in = and iM = 0, as if the system were autonomous 0. 
For this we assume the solutions: 

in = ci exp (At) 
Im = C2 exp (At) 

(3) 

iM = C3 exp (At) 
iE = C4 exp (At) 

drop the Heaviside 6*— functions by assuming cj > dj, and replace Q into 
equation (j2)). The characteristic equation associated to system Q is then 
obtained: 



(2) 
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<?rM (l - 



Pi {ci-di sin($)) 

(4) 

If all the roots of equation (jl)) have negative real parts, then the equilib- 
rium without disease is stable, that is, the origin is an atractor.. As shown 
in [4 , the first root that crosses the imaginary axis do so through the real 
axis and this happens when 



T3(f\ _ 2 Nmjt-Tj) aexp(-/JJ^fr)fec Snit) 

(5) 



+ 

pj{cj-dj sm<I>)grAf 

MA/(MB+P/(ci-rf/ sin*)) 



> 1 



Note that the first term in equation is exactly the expression proposed 
in |S| for the so-called 'basic reproduction number'. 

The intuition behind the above procedure is the following. System has 
'no-mass', that is, it responds to perturbations instantaneously. Therefore, 
we can find the time t at which the stability of the trivial solution of system 
(j2I), that is, iE = 0, in = and iM = becomes unstable. We have numeri- 
cally checked that the time t at which the trivial solution (no-disease) of the 
autonomous system becomes unstable {R > 1) corresponds approximately 
to the moment at which the epidemic takes off, that is, when the epidemic 
in system (P) begins to increase as a result of the introduction of a small 
amount of disease at time t = 0. 
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3 Qualitative analysis of the system's behav- 
iors 



In this section we analyze qualitatively all the possible behaviors of the sys- 
tem when a small amount of disease is introduced into a previously uninfected 
population and when R{t) is in its minimum value (winter time), that is, we 
ser = TT. We do so by using R(t) as in equation (0) and by numerically 
simulating system Q with parameters values as in table 1. 

Table 1 

The initial conditions were obtained by using the values of the carrying 
capacities (see table 1) and running the system without disease and taking 
the lowest values corresponding to the peak of the winter. The values are 
SniO) = 200,000, 5^/(0) = 85,000 and SsiO) = 930,000. The disease 
was introduced through a single infected egg, that is /^(O) = 1 and all the 
remaining variables equal to zero. 

We analyze two epidemiological scenarios, one in which R{t), in the ab- 
sence of infection, is most of the time above one, and another in which R{t) 
is most of the time below one. In the first case, if a small amount of infection 
is introduced we observe a pattern shown in figure 1. In the second case a 
small amount of infection introduced generates a pattern shown in figure 2. 

Figure 1 
Figure 2 

In figure 1 the intensity of transmission is relatively low (a = 3.7 days^^) 
and we see a first peak followed by a succession of outbreaks forming a 
damped oscillation pattern and the disease disappears. In other words, after 
the first outbreak the infection transmission decreases to levels inferior to 
that of the previous cycle. As the system oscillates the time interval during 
which R{t) > 1, that is, the system is above the threshold for transmission, 
is insufficient for keeping transmission, we have the pattern observed. 

In figure 2 the intensity of transmission is higher (a = 4.3 days~^) than 
that shown in figure 1 and, consequently the time interval during which 
R{t) > 1, that is, the system is above the threshold, is larger. In this case 
the amplitude of the consecutive outbreaks increases until the fraction of 
immune individuals reaches a herd immunity threshold and the disease dies 
out in a damped oscillation pattern. 
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We have numerically found that there is an increase in the amplitude of 
consecutive outbreaks with a = 4.3 days^^, whenever the preceding period 
of time R{t) is above the threshold for transmission is greater than a certain 
time interval (about 190 days). Indeed, with a = 3.7 days~^, as mentioned 
above, there is a first outbreak and the subsequent peaks formed a damped 
oscilation and the time interval R{t) is above the threshold for transmission 
is greater than 190 days only for the first peak. 

When we simulate the system with parameters that make R(t) < 1 for all 
times, the disease cannot invade the population and disappears exponentially. 

Those are the only three possible qualitative patterns generated by a 
small amount of disease introduced into an entirely susceptible population 
when R(t) is at its minimum value. Let us concentrate in the first two pat- 
terns, which can better be visualized in figures 3 and 4, where the threshold 
parameter R{t) is shown as function of time for each of the above cases. 

Figure 3 
Figure 4 

In figure 4 we can note the herd immunity effect acting after the second 
peak. 

Other interesting results are shown in figures 5 and 6, in which the time 
oscillation of the 'total amount of disease', d{t), defined as 

d{t) = ^{lH{t)f + (LM(t))' + {lM{t)f + {lE{t)f (6) 

is plotted together with R{t), for both cases of low and high intensities of 
transmission. It can be noted from the figure, that the points in which R{t) 
crosses 1 corresponds, approximately, to maximums and minimums of the 
function d{t). Note also that, in both cases the peaks and troughs of d{t) 
occur slightly after R{t) crosses 1, decreasing and increasing, respectively, as 
if the system has a small 'inertia'. 

Figure 5 
Figure 6 

4 Sensitivity analysis 

In this section we describe the sensitivity of the patterns to the amount of 
disease introduced at t = 0, and the sensitivity of the patterns to the time 
of the year at which the disease is introduced. 
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4.1 Sensitivity of the patterns to the amount of disease 
introduced 

It can be numerically cheked that when the transmission is relatively low, 
that is, when R{t) > 1 for periods of time less or equal to around 190 days, 
an increase in the amount of the disease introduced at t = changes the 
pattern shown in figure 1 by increasing the peaks almost linearly. In figure 
1 we introduced at t = one infected egg. When we introduced 5 infected 
eggs, the peaks amplitudes are multiplied by 5. Naturally, if a large amount 
of disease is introduced then herd immunity can distort the pattern. 

In the case when transmission is relatively high, that is, when R(t) is most 
of the time above one, an increase in the amount of disease introduced at 
t = can distort the pattern due to herd immunity. If the amount of disease 
introduced is sufficiently large we can observe just a single peak, that is, the 
disease disapear by a substantial decrease in the number of susceptibles. 

4.2 Sensitivity of the patterns to the time at which the 
disease is introduced 

We can vary the time of the year at which the disease is introduced by varying 
(f). If the disease is introduced when -R(O) < 1, then the disease dies out until 
the moment when R (t) crosses one from below. The rest of the development 
is identical to the pattern shown in figure 1 or 2, depending on the intensity of 
transmission. On the other hand, if the disease is introduced when -R(O) > 1, 
then the disease imediately increases and the subsequent pattern is as shown 
in figure 1 or 2, depending on the intensity of transmission. 

5 Final comments 

This paper presents a novel, as far as we know, approach to analyze the 
response of a non-autonomous system to a perturbation. This is quantified 
by an approximate expression for the threshold condition that determines 
whether the system will amplify or reduce a small of disease introduced at 
time t. We should warn the reader that we have no mathematical proof of the 
correctness of the threshold expression here deduced but intuition and the 
numerical investigation presented above suggests that it is basically correct. 
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A possible application exemplified in this paper is the case of dengue fever, 
in which a seasonal variation in the density of vector mosquitoes determines 
the intensity of transmission. In another paper we used a similar model 
to explain the question of overwintering, that is, how dengue fever survives 
through the winter's dry and cold season. 

Finally, we think that the analysis proposed in this paper could be applied 
to other vector-borne infections and also to some directly transmitted diseases 
that show seasonality in the intensity of transmission. 
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Captions for the Figures 



Figure 1. Number of infected humans for the case of low intensity of 
transmission (a = 3.7 days"^). There is a first outbreak resulting from the 
introduction of a small amount of infection in a previously uninfected popu- 
lation followed by a pattern of damped oscillation until the disease disappear. 
Simulation begins at the peak of the winter. 

Figure 2. Number of infected humans for the case of high intensity 
of transmission (a = 4.3 days~^). After the first outbreak resulting from 
the introduction of a small amount of infection in a previously uninfected 
population there arc subsequent outbreaks with larger amplitudes until herd 
immunity is achieved and the disease gradually disapears. Simulation begins 
also at the peak of the winter. 

Figure 3. The threshold R{t) in the case of low transmission (a = 3.7 
days~^). 

Figure 4. The threshold R{t) in the case of high transmission (a = 4.3 

days^-*^). 

Figure 5. The threshold R(t) and 'total amount of disease' d{t) in the 
case of low transmission (a = 3.7 days~^). The peaks and troughs of d{t) 
occur slightly after R{t) crosses 1. 

Figure 6. The threshold R{t) and 'total amount of disease' d{t) in the 
case of high transmission (a = 4.3 days~^). Again, the peaks and troughs of 
d{t) occur slightly after R{t) crosses 1. 



12 



Table 1 



Pcirameter 


Meaning 


Value 


a 


Average Daily biting rate 


see text 


h 


Susceptibility to Inlection 


0.1 




Humans Natural Mortality Rate 


4 X 10~^days~^ 




Humans iVlaltnusian Parameter 


1 days 


kH 


Humans Carrying Capacity 


10^ 




Dengue Induced Mortality m Humans 


10 ''days ^ 


1h 


Humans Recovery Rate 


0.143 days~^ 


Ps 


Proportion of non-infected eggs that reach adult phase 


0.15 


cs 


CUmatic factor modulating winters and summers 


0.08 


ds 


CUmatic factor modulating winters and summers 


0.06 


f 


Frequency of the seasonal cycles 


2.8 X 10"^days" 




Natural mortality rate ot mosquitoes 


0.263 days ^ 


T 


Extrinsic incubation period of dengue 


7 days 




Dengue induced mortality in mosquitoes 


negligible 


tm 


Eggs Malthusian parameter 


50 days~^ 


Pi 


Proportion of infected eggs that reach adult phase 


0.15 


Cl 


Climatic factor modulating winters and summers 


0.06 


di 


Climatic factor modulating winters and summers 


0.06 


9 


Proportion of infected eggs laid by infected females 


see text 


kE 


Eggs Carrying Capacity 


10*^ 




Natural mortality rate of eggs 


0.1 days-i 
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